Low-complexity quantum codes designed via codeword-stabilized framework 
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We consider design of the quantum stabilizer codes via a two-step, low-complexity approach based 
on the framework of codeword-stabilized (CWS) codes. In this framework, each quantum CWS code 
can be specified by a graph and a binary code. For codes that can be obtained from a given graph, 
we give several upper bounds on the distance of a generic (additive or non-additive) CWS code, 
and the lower Gilbert- Varshamov bound for the existence of additive CWS codes. We also consider 
additive cyclic CWS codes and show that these codes correspond to a previously unexplored class 
of single-generator cyclic stabilizer codes. We present several families of simple stabilizer codes with 
relatively good parameters. 



I. INTRODUCTION 



It was the invention of quantum error correcting codes 
lUHl (QECCs) that opened quantum computing (QC) as 
a theoretical possibility. However, high precision required 
for error correction [^]-|l3| combined with the large num- 
ber of auxiliary qubits necessary to implement it, have 
so far inhibited any practical realization beyond proof- 
of-the-principle demonstrations Jl4|-p0| . 

In any QECC, one needs to perform certain many-body 
quantum measurements in order to decide how to cor- 
rect the encoded state. The practical difficulty is that a 
generic code requires measurements which are both com- 
plicated and frequent at the same time. It is therefore 
clear that a quantum computer can only be build via a 
thorough optimization at every step of the design. In 
particular, code optimization targets codes that combine 
good parameters with fairly simple measurements. It is 
also desirable to parallelize these measurements given a 
specific on-chip layout of a QC architecture. 

To date, the main focus of the QECC-research has been 
on finding good codes with the traditional code parame- 
ters, which are the block length n, code dimension K, and 
code distance d (or code rate R = (log 2 K) jn and the rel- 
ative distance 5 = d/n). For stabilizer codes ^2|, we 
also consider the number of encoded qubits k = log 2 K. 

A number of stabilizer codes |23| have been designed 
that meet or nearly achieve the existing bounds on dis- 
tance d for the given k and n. Code parameters can be 
further refined by going beyond the family of stabilizer 
codes. One example is a recently introduced framework 
of codeword-stabilized (CWS) quantum codes p4j-g7|. A 
qubit CWS code Q = (Q,C) (in standard formjis deter- 
mined by a graph Q and a classical binary code C. CWS 
codes include all stabilizer codes as a subclass (the cor- 
responding binary code C must be linear), but also the 
codes which have been proved to have parameters supe- 
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rior to those of any stabilizer code] 
nately, typical gains in code dimension K correspond to 
a fraction of a qubit. Moreover, error-correcting algo- 
rithms known for general non-additive CWS codes have 
exponential complexity ]3l| p4[| , as opposed to polynomial 



complexity of the stabilizer codes. 

Even for the relatively simple additive codes, their op- 
timization is a very difficult problem that has exponen- 
tial complexity. This is one of the main reasons as to 
why the two relatively simple code families are almost 
exclusively used among stabilizer codes to estimate the 
threshold accuracy required for scalable quantum compu- 
tation: the concatenated codes [Q, ||, |7|-|l2| and the surface 
codes H |l3| which originated from the toric codes [|35|. 
Both families have very low code rates that scale as in- 
verse powers of code distance. 

In this work we explore how the framework of CWS 
codes can be used to relegate the design of quantum sta- 
bilizer codes to classical binary linear codes in order to 
simplify the overall design. In particular, we formulate 
several theorems framing the parameters of an additive 
CWS code which can be obtained from a given graph. 
We also suggest a simple decomposition of the F4 gen- 
erator matrix corresponding to the stabilizer in terms of 
the graph adjacency matrix and the parity check ma- 
trix of the binary code. Finally, we design several graph 
families corresponding to regular lattices which result in 
some particularly good codes. These include graphs with 
circulant adjacency matrix which can be used to con- 
struct single-generator cyclic additive codes, a class of 
codes overlooked in previous publications. In particu- 
lar, we prove the existence of single-generator cyclic ad- 
ditive codes with the parameters [[km, k, m]], k > 10 
and [[t 2 + (t + l) 2 ,l,2t + 1]] (version of toric codes). 
Note that these code families have distances that are 
not bounded, unlike any CWS code families constructed 
previously p5|, |3f| . 

The paper is organized as follows. In Sec. [H], we intro- 
duce the notations and briefly review some known results 
for quantum and classical codes. In Sec. pi] , we establish 
several upper bounds on general CWS codes. In Sec. |^ 
we give a CWS decomposition of the F4 matrix corre- 
sponding to the stabilizer generators. In Sec. [y], we for- 
mulate the Gilbert- Varshamov (GV) bounds for additive 
CWS codes which can be obtained from a given graph. 
Cyclic additive CWS and more general single-generator 
additive cyclic codes are considered in Sec. VI where we 
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discuss their properties an d giv e several examples. We 
give our conclusions in Sec. VI]. 



form the n-qubit Pauli group & n of size 2 2n+2 , 



II. NOTATIONS AND SOME KNOWN RESULTS 
A. Classical and quantum error correcting codes 

A classical g-ary block error-correcting code (n, K, d) q 
is a set of K lcngth-77, strings over an alphabet with q 
symbols. Different strings represent K distinct messages 
which can be transmitted. The (Hamming) distance be- 
tween two strings is the number of positions where they 
differ. Distance d of the code C is the minimum distance 
between any two different strings from C. 

In the case of linear codes, the elements of the al- 
phabet must form a Galois field F q ; all strings form n- 
dimensional vector space F™. A linear error-correcting 
code [n,k,d] q is a fc-dimensional subspace of F™. The 
distance of a linear code is just the minimum weight of 
a non-zero vector in the code, where weight wgt(c) of a 
vector c is the number of non-zero elements. A basis of 
the code is formed by the rows of its generator matrix G. 
All vectors that are orthogonal to the code form the cor- 
responding (n — fc)-dimensional dual code, its generator 
matrix is the parity-check matrix H of the original code. 

For a binary code C[n, k, d], the field is just F 2 = {0, 1}. 
For a quaternary code C, the field is F4 = {0, l,w,ZZ?}, 
with 



uT = uj + 1, 



; 3 = 1, and uj = uj 2 . 



(1) 



For non-binary codes, there is also a distinct class of 
additive classical codes, defined as subsets of F™ closed 
under addition (in the binary case these are just linear 
codes). A code C is cyclic if inclusion (cq, ci, . . . , c„_i) G 
C implies that (c„_i, Cq, ci, . . . , Cn—2) S C. Codes that 
are both linear and cyclic are particularly simple: by 
mapping vectors to polynomials in the natural way, c — > 
c(x) = cq+c\x+. . .+c n -\x n ~ , it is possible to show that 
any such code consists of polynomials which are multiples 
of a single generator polynomial g(x), which must divide 
x n — 1 (using the algebra corresponding to the field ¥ q ). 
The quotient defines the check polynomial h(x), 



h{x)g(x) = x n - 1 



(2) 



which is the generator polynomial of the dual code. The 
degree of the generator polynomial is dcgg(x) = n — 
k. The corresponding generator matrix G can be chosen 
as (the first k rows of) the circulant matrix formed by 
subsequent shifts of the vector that corresponds to g (x). 

Qubit quantum error correcting codes are defined in 
the complex Hilbert space T-Lf n , where H2 is the Hilbert 
space of a single two-level system. % 2 is formed by all 
vectors a |0) + (3 |1) with a,/3 G C, and the inner prod- 
uct such that the two states are orthonormal, = 
i, j G {0, 1}. Any operator acting in Hf n can be repre- 
sented as a linear combination of Pauli operators which 



1 {I,X,Y,Z} 



0, ...,3 



(3) 



where A, Y, and Z are the usual Pauli matrices, and 
I is the identity matrix. The weight wgt(-E) of a Pauli 
operator E is the number of non-identity terms in the 
corresponding tensor product. 

All Pauli operators are unitary; they are also Hermi- 
tian with eigenvalues ±1 when the phase factor i m in 
Eq. (||) is real-valued, m = 0, 2. A state \tp) is stabilized 
by a Hermitian Pauli operator M if M = \ip). A lin- 
ear space Q is stabilized by a set of operators M. if each 
vector in Q is stabilized by every operator in M.. 

An ((n,K,d)) quantum error- correcting code is a K- 
dimensional subspace of the Hilbert space T-lf n . Such 
a subspace can be described by an orthonormal basis 
{K)}|=i- Let £ C be some set of Pauli errors. A 
QECC detects all errors E G £ if and only if M, 



(j \E\ i) = C E 8 l3 



(4) 



where Ce only depends on the error E, but is indepen- 
dent of the basis vectors. A QECC has distance d if it can 
detect all Pauli errors of weight (0? — 1), but not all errors 
of weight d. The errors in the set £ can be corrected if 
and only if all the nontrivialpairwise combinations of er- 
rors from £ are detectable |j| . Thus a distance-G? code 
corrects all errors of weight s <t= [(d — 1)/2J. 

The code Q is non-degenerate if linearly-independent 
errors from £ produce corrupted spaces E(Q) = {E : 
G Q} that are linearly independent. Otherwise, the 
code is degenerate, implying the existence of at least 
two mutually degenerate linearly independent operators 
{£^1,^2} G £ which act identically on Q. 

The code is called pure if linearly independent errors 
from £ produce corrupted spaces that are not only lin- 
early independent, but also mutually orthogonal. For all 
codes considered in this work, non-degenerate codes are 
also purep^, |3^| . 

Two codes are considered equivalent if they differ just 
by qubit order, and/or discrete rotations leaving each of 
the singlc-qubit Pauli groups invariant. The latter are 
called local Clifford (LC) transformations. 



B. Stabilizer quantum error correcting codes 

Here we briefly review the well-known family of sta- 
bilizer codes pij . An [[n, k,d]] stabilizer code Q is a 2 k - 
dimensional subspace of the Hilbert space Hf n stabilized 
by an Abclian group C & n with n — k Hermitian Pauli 
generators, = (Gi, . . . , G n -u). Explicitly, 



Q = {\ip) :S\ip) = \ip),vs g y} 



(5) 



Such a code exists only if —1 ^ 5?. The group 5? is 
called the stabilizer of the code. Changing the sign(s) of 
one or several of the generators Gi results in replacing Q 
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with one of 2 n ~ k — 1 equivalent codes whose direct sum 
(together with Q) is the entire space @> n . 

The normalizer of 5? is a set of Pauli operators gener- 
ating unitary transformations that leave y invariant, 



Jf = {U e &> n : U^SU = S,VS e y} 



(6) 



Elements of the normalizer form a group commuting with 
5? but not necessarily with each other. It is possible 
to construct 2k logical operators Xj, Zj, j = l,...,k 
belonging to with the usual commutation relations 
that generate the normalizer when the generators of 
are included §^J§. The Abelian sub group of jY , S^q = 
(Gi, . . . , G n -k, Zi, . . . , Zk), becomes a maximal Abelian 
subgroup of £P n when the generator il is also included. 

The group o5^o stabilizes a unique stabilizer state | s) = 
1 . . . 0) , an [[n, 0, d']] stabilizer code, while the operators 
Xj generate the basis of the code, i.e., 



ci 



■ Cfe) 



Xi 



■ K k \s) 



(7) 



By convention, the stabilizer state is considered non- 
degenerate, and its distance d! is the minimum weight 
of a non-trivial member of the group J^ . 

For stabilizer codes, phases of (Hermitian) Pauli op- 
erators are only needed to choose one of the equivalent 
codes in Eq. (^|) , as well as to introduce the commutation 
relations. It is convenient to drop the phases and map 
the Pauli operators to two binary strings, v, u £ {0, 1}™ 

il, 



u 



X V Z U -> (v, u) 



(8) 



where X v = X^X% 2 . . . X^ , Z u = Z^ 1 Z^ . . . Z^ , and 
m! = 0, . . . , 3 is generally different from that in Eq. (J3j) . 
This map preserves the operator algebra, with a product 
of two Pauli operators U\ and U2 corresponding to a sum 
of the corresponding binary vectors (vi, ui) and (V2, U2). 

The map (^) can be taken one step further Q to qua- 
ternary codes, by introducing F4 vectors e = u + wv [see 
Eq. (|l|) ; note that this mapping differs slightly from that 
in Rcf . p2| . We will denote this combined map as a func- 
tion e = 4>(U). Note that up to a phase this association 
also allows us to define <j>~ 1 (e). To be specific, for the 
Pauli operator 4>~ 1 (e) we will set ml = v • u in Eq. (JsJ) , 
which corresponds to m — in Eq. (||). 

It is easy to check that two Pauli operators commute 
if and only if the symplectic scalar product vi • 112 + 
Ui • V2 vanishes (mod 2). In terms of the corresponding 
{ei,e2} C F4, this corresponds to the vanishing of the 
trace inner product 



ei * e 2 = ei • e 2 + ei • e 2 , 



(9) 



where ej = Ui + uJv.;, i = 0, 1. 

A dual code to an additive F4 code C (equipped with 
trace inner product) is defined as [^| 

C± = {e' £FJ:e'*e=0, for all e G C} . (10) 



If C C C±, one says C is self-orthogonal. A classical 
additive code C corresponding to a set of operators 5?\ 
is self-orthogonal if and only if oS^i is an Abelian group. 
Thus any quantum stabilizer code can be described as a 
self-orthogonal classical additive code over F4. The fol- 
lowing theorem is applicable to additive F4 codes (variant 
of Theorem 2 from Ref. @): 

Theorem 1. Suppose C is an additive self-orthogonal 
code in FJ, containing 2 n ~ k vectors, such that there are 
no vectors of weight < d in C± \ C . Then </> _1 (C) de- 
fines a stabilizer of an additive QECC with parameters 
[[n,k,d]]. 

Example 1. The well-known Calderbank-Shor-Steane 
(CSS) [[7,1,3]] code [7^/ has the stabilizer with the 
generators]^! 

XXXXIII, XXIIXXI,XIXIXIX, 
ZZZZIII, ZZIIZZI, ZIZIZIZ, (11) 

and the logical operators 



X = zzzzzzz, z = xxxxxxx. 



(12) 



As any CSS code, this code is linear. Qubit permuta- 
tions also give an equivalent cyclic linear code with the 
generator polynomial g(x) = 1 + x + x 2 + x 4 ; g(x) is a 
factor ofx 7 -l. The corresponding check polynomial is 
h(x) = 1 + x + x 3 . 



C. Codeword stabilized codes 

CWS codes |2{| represent a general class of nonaddi- 
tive QECCs. A general CWS code is defined in terms of 
a stabilizer state \s) and a set of K mutually commut- 



ing codeword operators W 
[Cf. Eq. @], 



W}<=i C 9> n . Explicitly 



Q = span({W i \s)}g 1 ). 



(13) 



For non-trivial CWS codes, this construction coincides 
with union-stabilizer (USt) codes [[II], restricted to the 
zero-dimensional originating code. 

Any stabilizer state is LC-equivalent to a graph state 
p2| E5[ , a stabilizer state with the stabilizer group Sfg = 
(Si, . . . ,S n ) whose generators Si are determined by the 
adjacency matrix R G {0, l}" xn of a (simple) graph Q, 



S, = \,/ v 



(14) 



where r^, i = 1, . . . , n denotes the i-th row of R. In fact, 
such a graph is usually not unique, even after account- 
ing for graph isomorphisms. The full set of LC-equivalent 
graph states can be generated by a sequence of local com- 
plementations, operations on a graph where the subgraph 
corresponding to a neighborhood of a particular vertex 
is inverted. Such graphs are called locally equivalent [^6] . 

Any CWS code ((n,K,d)) is LC equivalent to a CWS 
code in standard form, defined by an order-n graph Q 
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and a classical binary code C containing K binary words. 
The graph defines the graph state, while the vectors of 
the classical code c, € C are used to generate the code 
word operators, W, = Z c ' . Thus, 

Q = S pM{Z c '\ S )}f =1 ). (15) 

It is customary to use notation Q = (Q, C) for CWS codes 
in standard form. 

The key simplification of the CWS construction comes 
from the fact that the basis states Wi \s) are eigenvectors 
of the graph stabilizer generators, 

SiWi \s) = ±W t \a) , 5, 6 S>g. (16) 

Thus, a Pauli operator in the form (|J) can be transformed 
to a Z-only operator Z Glg ^ u \ where the graph image of 
the operator U is the binary vector 

n 

C\g(U) Eu + ^v i r 1 (mod2). (17) 

i=l 

The error correcting properties of a quantum CWS 
code Q = (G,C) and the classical code C are related by 
the following 

Theorem 2. (after Theorem 3 from Ref. Jgjj/j Con- 
sider a CWS code Q = (G,C). An error E such that 
Clg(-E) 7^ 0, is detectable in Q if and only if the binary 
vector Clg(E) is detectable within the code C. An error 
E such that Gig (E) = is detectable in Q if and only if 
Z C E = EZ C for allceC. 

The case Clg(E) ^ corresponds to pure (non- 
degenerate) errors, while C\g(E) = indicates that the 
error is in the graph stabilizer group S*g ; the correspond- 
ing dctcctability condition is a requirement that the error 
must be degenerate. 

While in general CWS codes are non-additive, they 
include all stabilizer codes as a subclass. A CWS code 
Q = (G,C) is additive if C is a linear code[^5|. The stabi- 
lizer 5? of an additive CWS code in standard form is 
a subgroup of the graph stabilizer 5?g\ it can be ob- 
tained from the graph-stabilizer generators by a sym- 
plectic Gram-Schmidt orthogonalization procedure [|34|. 
Conversely, the representation (Q) of an additive code 
corresponds to a general CWS code; an LC transforma- 
tion may be needed to obtain the corresponding standard 
form, and one can always find a standard form where C 
is linear In the following we will always assume such 
a representation. 

Example 2. The smallest single- error- correcting code is 
the linear cyclic [[5, 1,3]] corfe/y, \^\ , [J^/ with the gener- 
ator polynomial g{x) = 1 + ujx + ujx 2 + x 3 which divides 
x — 1. This code is unigue; its stabilizer generators can 
be obtained as cyclic permutations of a single operator, 
XZZXI , and the logical operators are 

X = ZZZZZ,~Z = XXXXX. (18) 




FIG. 1. Left: 5-ring graph corresponding to the [[5, 1, 3]] code 
in Example Ex. |j| Center: "Kite" graph corresponding to the 
degenerate code [[6, 1,3]] from Example [j. Right: The graph 
corresponding to the cyclic [[7, 1,3]] code from Example^. 

The corresponding CWS codejfefy can be generated from 
the 5-ring graph in Fig. |t| (left), and the binary code has 
a single generator c = (11111). Note that both the graph 
and the binary code preserve the original cyclic symme- 
try. 

Example 3. There exist only two ineguivalent single- 
correcting codes [[6,1,3]]; both are degenerated^ . One 
of the codes is obtained from the code in Example by 
adding a gubit; the graph of the corresponding CWS code 
can be chosen as a 5-ring [Fig. ^| left] and a discon- 
nected vertex i = 6; the binary code C is generated by 
c = (111110). The degeneracy group is generated by 
Sq = Xq. The stabilizer generators for the second code 
are listed in Ref. J#^/. This code corresponds to the graph 
in Fig. U (center), while the binary code is generated by 
Ci = (011100). While there are three bits which are not 
involved with the classical code, they cannot be dropped 
as they are part of the entangled state. The degeneracy 
group is generated by S1S2 = X1X2 (the equivalence fol- 
lows from the fact that the first two vertices of Q share 
all of their neighbors). 

Example 4. The linear cyclic [[7,1,3]] CSS code from 
the Example M is LC equivalent to a CWS code with the 
graph in Fig.\^ (Right). The corresponding classical code 
is given by C = {0000000,1110000}. Note that neither 
the graph nor the binary code is explicitly symmetric with 
respect to cyclic permutations of the qubits. Note also 
that an inequivalent CWS cyclic [[7, 1,3]] code exists; we 
constructed such a code among others in Example |7J, see 
Table @. 

III. UPPER BOUNDS FOR CWS CODES 

In this section we give upper bounds on general CWS 
codes in terms of the properties of the corresponding 
graph G and the binary code C. 

Lemma 1. The distance d of the CWS code Q = (G,C) 
cannot exceed that of C . 

Proof. Indeed, any "classical" error in the form E = 
i m Z u is mapped by Eq. ( p"7| ) to the binary vector u. If 
E is detectable by Q, u should be detectable by C. □ 
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Lemma |l| concerns with errors which are dealt with by 
the binary code. On the other hand, a CWS code is an 
enlargement of the code formed by the graph state. The 
following observation has been made in Ref . ||^] : 

Lemma 2. The distance d of a nondegenerate CWS code 
((n,K,d)) is limited by the distance d'(G) of the graph 
stabilizer state, d < d'(Q). 

It follows from the fact that any member of the graph 
stabilizer is either a degenerate error, or it is a non- 
detectable error. Note, however, as illustrated by the 
Example [| in general, the distance of a CWS code can 
actually be bigger than that of the graph stabilizer state. 

For a binary code C, we will say that the j-th bit is 
involved in the code if there are vectors in the code for 
which the value of j-th bit differ, c\ ^ c? 2 . Alternatively, 
if the all-zero vector is in the code (which can always 
be arranged), the condition is that there is a vector c G C 
where j-th bit is non-zero, c J ^ 0. 

Lemma 3. For a CWS code Q = (Q,C) with K > 1, let 
us assume that j-th bit is involved in the code C. Then the 
graph- stabilizer generator Sj violates the error detection 
condition in Theorem @. 

Proof. Since the generator Sj is in the graph stabilizer, 
Clg(Sj) = 0, one has to check the degenerate condition 
in Theorem |2|. The commutativity of Sj with a given 
Z c is determined by the j-th bit of c; conditions of the 
Lemma ensure that only one of the two vectors commute 
with Sj. □ 

Note that this means that the code distance cannot 
exceed that of any Sj corresponding to bits involved in 
the binary code. Since at least d bits must be involved 
in the binary code, Lemma [j] guarantees the following 
bound 

Theorem 3. The distance d of a CWS code Q = (Q, C) 
cannot exceed the d-th largest weight of Si, minimized 
over all graphs that are locally- equivalent to Q. 

We will also be using the following 

Corollary 3.1. For a graph Q with all vertices of the 
same degree r, the distance of a CWS code Q = (Q,C) 
cannot exceed r + 1 . 

In particular, for any ring graph r = 2, which gives 
d < 3, for any double-ring graph|26| d < 4, and for a 
large enough square lattice wrapped into a torus, d < 5. 

Obviously, to maximize the distance of a CWS code, 
one may want to maximize the distance of the binary 
code C. To this end, it is a good idea to make sure that 
every bit is involved in C. For such codes, we have 

Theorem 4. The distance of a CWS code Q = (Q, C) 
where the binary code C involves all bits cannot exceed 
the minimum weight of Si, minimized over all graphs that 
are locally-equivalent to Q . 



IV. ADDITIVE CWS CODES AND 
QUATERNARY CODES 

The stabilizer of an additive CWS code Q = (Q, C) is 
a subgroup of the Abclian graph stabilizer S^q, and its 
generators Gi G 5? can be expressed as products of graph 
stabilizer generators Si Explicitly, 

n 

G t = l[Sp\ (19) 

3=1 

where P G {0,l}"~ fcxn is the corresponding matrix of 
binary coefficients. With the help of Eq. (|l4]), we obtain 
the following decomposition for the generator matrix G 
of the associated additive F4 code C, 

G = P{ul + R), (20) 

where R G {0,1}" XTI is the symmetric graph adjacency 
matrix with zeros along the diagonal, and 1 is the n x n 
identity matrix. The relation between the binary code 
C and the quaternary code C is explicitly given by the 
following 

Lemma 4. The additive F4 code C with the generator 
matrix is the map 4 > {^') °f ^ e stabilizer 5? of the 
additive CWS code Q = (Q, C) generated by the graph 
with the adjacency matrix R and the linear binary code 
C if and only if P is the parity check matrix of C. 

Proof. Use the basis vectors in the form ( p"5| ) and the 
commuting operators ( |l9|) corresponding to the rows of 
the matrix G [Eq. (p0|)] . Direct calculation gives 

n 

Gi Z c \s) = Y[ Sf' J Z c \s) = (-l) p ^ cJ Z c \s) , (21) 
i=i 

were summation over repeated indices is assumed. The 
statement of the Lemma (both ways) follows from the 
definition of the parity check matrix. □ 

We can now easily relate the error detection conditions 
for additive codes in Theorems [l] (codes over F4) and || 
(CWS codes). The code C in Theorem |l| is given by 
additive combinations of the n — k rows of the generator 
matrix (p0|). Evaluating the outer trace product of the 
generator matrix (^o|) with a vector e = u+wv, we obtain 
the condition for the vector to be in Cj_ 

= G*e = P(u + Rv). (22) 

This uniform linear system of n — k equations with 2n 
variables has n + k linearly-independent solutions. The 
corresponding basis can be chosen as a set of k linearly- 
independent "classical" vectors e^ = u^, where Pui = 
and the corresponding Vi = 0, i = 1, . . . , k, plus n 
linearly-independent vectors such that Uj = Rvj, j = 
k + 1, . . . , k + n. Some linear combinations of the latter 
vectors are actually in C. These can be found using the 
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identity (Cj_)j_ = C: the corresponding Vj have to satisfy 
u, • Vj = 0, which precisely corresponds to the degenerate 
case, C\g(E) = 0, in Theorem 0. 

General theory of CWS codes guarantees that gener- 
ator matrix of a quantum code equivalent to any addi- 
tive self-orthogonal code over F4 can be decomposed in 
the form (ppf). Conversely, any matrix in the form ( |20| ) 
with binary matrices P and R generates a self-orthogonal 
code over F4 as long as the matrix R is symmetric. We 
use it in the following section to prove the lower Gilbert- 
Varshamov (GV) bound for the parameters of an additive 
CWS code which can be obtained from a given graph. 



V. GV BOUND FOR THE ADDITIVE CWS 
CODES WITH A GIVEN GRAPH 

The GV bound is a counting argument which non- 
constructively proves the existence of codes with param- 
eters exceeding certain threshold. The argument is based 
on the fact that the set of possible codes (vector spaces) 
vastly outnumbers the set of vectors. Then, if we count 
all codes of a given length n, and then subtract the num- 
ber of codes that contain any vector of weight d — 1 or 
less, the remaining codes (if any) will all have distance 
d or more. This "greedy" argument ignores any possi- 
ble double counting of codes that contain several small- 
weight vectors. Note that the GV bound necessarily gives 
asymptotically good codes with relative distance d = d/n 
and code rate R = k/n bounded away from as n — > 00. 

For the entire class of pure stabilizer codes, the asymp- 
totic GV bound |49f| states that there exist codes of length 
n — >• 00 such that 



S\og 2 3 + H 2 (6) >1-R, 



(23) 



where H^S) = — <51og 2 5— (1 — 5) log 2 (l — 5) is the binary 
entropy function. We are going to prove that the same 
bound also holds for pure CWS codes corresponding to 
a given graph Q, as long as d < d'{Q) [see Lemma g]. We 
are using Eq. ( |20| ) to parameterize the stabilizer matrices; 
the resulting F4 codes are automatically self-orthogonal. 
Let N^ k be the number of CWS codes that have length 
n, dimension at least k, and correspond to a given graph 
Q. Let also k be the number of such codes which 
contain a given vector e = u + wv, wgte < d (Q). in 
C± [sec Eq. (|2(i|)1. The corresponding condition ( |To| ) is 
given by the trace inner product (p2[). For wgte < d'(Q), 
the binary vector c = u + i?v is always non-zero (which 
also guarantees that e £ C). As a result, N^ nk and 

k represent the corresponding numbers for the binary 
codes. 

Then the standard counting arguments p0[| show that 



(2 n - = ( 2 " - 



(24) 



1 vectors c belongs 

to the same number N2„ u of binary codes; also each of 



Here we use the fact that each of 2 

e,n.k 



number of quaternary vectors of weight s is 3 s (™) . Thus, 
for any graph Q, there exists a distance-^ CWS code as 
long as 



d-x 



s=l 



(25) 



Now we see that there exist [[n, k, d]} CWS codes for the 
graph Q with distance 

d = mm{dGv, d max }, (26) 
where d max is the distance of the graph state d'(Q) and 



IGV 



d-1 

max d : 3 * 

s=l 



< 



2" - 1 
2 fc - 1 ' 



(27) 



Note that thus obtained quantum codes are always 
pure, since the summation in the l.h.s. can only be ex- 
tended up to e? max — 1. Apart from this latter condition, 
Eq. (B7 | ) i s identical to the quantum Gilbert- Varshamov 
boundpsl for pure stabilizer codes, and takes the asymp- 
totic form ( p3| ) as n — > 00. 

The exact GV bound d > day for pure stabilizer codes 
(without the restriction on the distance) is recovered if 
we go over different graphs. Indeed, the GV bound ( |23| ) 
also applies for the special case of k = 0, corresponding 
to stabilizer states or self-dual codes [Ell. The GV bound 
on the relative distance is monotonous in k and reaches 
its maximum at 



5 k=0 w 0.189. 



(28) 



Then for given n and S < 6k=o one can always find a 
suitable graph such that the GV bound d < dav becomes 
more restrictive than the condition d < d maK . 

In practice, graphs with large distance d'(Q) are com- 
plicated (have too many edges). It is much easier to come 
up with graph families corresponding to a fixed graph- 
state distance d'(Q). For n — > 00, the corresponding code 
families approach the maximum rate R = 1 and have 
asymptotic redundancy r = n — k defined by Eq. (27): 



r < dlog 2 3 + nH.2(d/n). 



(29) 



N^ k binary codes contains 2 k 



1 nonzero vectors c. The 



It is readily verified that the r.h.s of estimate (g9() has 
the order of dlog 2 (3n/d) if d=const and n — > 00. 

Example 5. Graphs in the form of sufficiently large 
finite square lattice fragments [Fig. ^ have maximum 
distance c? max = 5, but this requires that the bits in 
the corners and around the perimeter be not involved in 
the classical code. Somewhat better redundancy can be 
achieved by avoiding only the bits in the corners, which 
gives d max = 4. Consider the family of classical codes 
where the codewords are obtained by taking all trans- 
lations of the pattern shown in Fig. with open sym- 
bols. The weight of any linear combination of such code- 
words is at least 4- The lattice shown corresponds to the 
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FIG. 2. Square-lattice additive CWS code [[25,4,4]]. Circles 
represent the qubits. All k = 4 translations of the empty- 
circle pattern form the classical codewords Cj. The generators 
of the stabilizer are formed as the products of the graph- 
state stabilizer generators along the directions parallel to the 
dashed lines. 



code [[25,4,4]], while general m x x m y lattice gives the 
code with the parameters [[m x m y , (m x — 3)(m y — 3), 4]]. 
Asymptotically, the redundancy for m x = m y is n — k w 
671 1 ' 2 , n — > 00, which is bigger than the logarithm in 
Eq. $2f\). However, the fraction of auxiliary qubits van- 
ishes as oc 1/n 1 ' 2 for large n. The code distance can be 
increased with higher- dimensional generalizations, e.g., 
we can generalize this construction to D '-dimensional hy- 
percubic lattice with 2D nearest neighbors so that the dis- 
tance is d = 2D in full analogy with the two dimensional 
case. The corresponding redundancy will scale with the 
area of the boundary. 



While the code in Fig. || serves as a good illustration of 
the concept of lattice codes with simple stabilizer struc- 
ture, it is still far from optimal. On the 5x5 square 
lattice we have constructed numerically a code [[25, 9, 5]] 
with weight- 7 codewords which can be mapped into each 
other by translations and rotations. This design is only 
one logical qubit short of the best-known generic code 
[[25,10,5]]. 



Example 6. Consider graphs in the form of Lx L square 
lattices wrapped into tori due to periodic boundary con- 
ditions. For L > 5, these graphs have the distance 
d'{Q) = 5. GV bound (||) shows that the CWS codes 
with the following parameters can be obtained for these 
graphs: [[25,4,5]], [[36,13,5]], [[49,24,5]], [[64,38,5]], 
[[81,53,5]], .... 



Example 7. Consider graphs in the form of L x L tri- 
angular lattices wrapped into tori due to periodic bound- 
ary conditions. These graphs have the distance d'(G) = 
min(L, 7). GV bound $Zb\) shows that the CWS codes with 
the following parameters can be obtained for these graphs: 
[[36,9,6]], [[49,15,7]], [[64,28,7]], [[81, 43, 7]], . . . . 



VI. SINGLE-GENERATOR ADDITIVE CYCLIC 
CODES 

Example [| shows that a cyclic additive code does 
not necessarily preserve its symmetry when converted to 
CWS standard form. By a cyclic additive CWS code 
we just mean a code which is cyclic in standard form, 
with a circulant graph. For such a code, Eq. ( pfj|) can be 
rewritten as the generator polynomial, 



g(x) = p(x)[u + r(x)], 



(30) 



where the polynomials p(x) and r{x) are binary, p(x) is 
the parity-check polynomial of a binary cyclic code (and 
therefore must divide x n — 1), while r(x) corresponds to 
a symmetric circulant matrix, 



r(x n - 1 ) = r(x) (mod / - 1). 



(31) 



Any such symmetric polynomial r(x) leads to a valid self- 
orthogonal additive code. The dimension of the quantum 
code corresponding to the generator polynomial (|o|) is 
k = n — degp(x). 

Previously, the additive cyclic QECCs were introduced 
in Theorem 14 of Ref. |22], stating that any such code 
has two generators. A smgfe-generator additive code de- 
scribed by Eq. ( |30| ) represents a new setting, in which the 
second generator is equal to zero. This c ondit ion gives 
a self-orthogonal additive code C [see Sec. II B] with no 
binary codewords (any e S C . e = u + low, has v^O). 

A somewhat wider class of single- generator cyclic ad- 
ditive codes can be also defined via Eq. (|30"|), without 
requiring the symmetry ( |3~l| ) of r(x). Then two codes Q 
and Q! that have a generator polynomial in the form ( |30| ) 
with the same p(x) = p'(x) are equivalent if and only if 



r{x) = r' (x) mod q(x), 



(32) 



where q(x) = (x n — l)/p(x) is the generator polynomial 
of the binary codej2^]. Such a polynomial ( |30] ) generates 
a self-orthogonal F4 code if and only if p2fl 



p(x)p(x n l )r(x 



n-l\ 



p(x)p(x n 1 )r(x) (mod x r 



(33) 

This guarantees self-orthogonality for any r(x) as long as 



p{x)p(x n - 1 ) = mods? 



1. 



(34) 



An alternative formulation of this sufficient condition is 
that the corresponding generator polynomial q(x) must 
contain no more than one root from each pair (a,cv _1 ) 
of mutually conjugate nth roots of unity, a" = 1. In 
particular, a self-reciprocal (palindromic) polynomial p2j , 



(35) 



always contains roots in pairs a and a -1 . For such poly- 
nomials Eq. (HJ) always fails (including the special case 
of q(x) = 1 + x which has only one root a = oT x = 1). 
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A. Single-generator cyclic codes from a binary code 

The algebraic condition (|^) on check polynomials for 
linear cyclic codes makes them simpler to implement but 
also dramatically restricts their number. In particular, 
the general counting approach to finding CWS codes 
[see Sec. |y[], where one first chooses a graph, and then 
searches for a suitable binary code can hardly be ap- 
plied to cyclic CWS codes. Even for classical binary 
cyclic codes, there are no counting arguments known 
to date that yield asymptotically good codes, let alone 
the stronger GV bound (see Research Problem 9.2 in 
Ref. |5(]). Also, long BCH codes — one of the major sub- 
classes of cyclic codes — are asymptotically bad and have 
a slowly declining relative distance S ~ (2 In R" 1 )/ log 2 n 
for any code rate R. On the other hand, binary cyclic 
codes often achieve the best known parameters (exceed- 
ing the GV bound) on short lengths n < 256. Thus, 
using simple cyclic codes in quantum design can yield 
both good parameters and feasible implementation on 
the short blocks. 

To better evaluate code distance of single-generator 
quantum cyclic codes (|30|), we will modify our count- 
ing approach of Sec. |Vj and begin with a binary cyclic 
code. Namely, we will fix some parity-check polynomial 
p(x) with a desired degree k among the binary factors 
of x n — 1. Then we will search for a polynomial r(x), 
either corresponding to a cyclic graph [see Eq. (|3l"|)], or 
satisfying the more general orthogonality condition (p3|). 
However, this transition will show that the parameters 
of quantum codes generated this way strongly depend on 
the chosen binary code. We will concentrate exclusively 
on the binary codes with irreducible generator polyno- 
mial q(x). We will show that the distance of such a cyclic 
CWS code is limited from below by the GV bound (or 
the variants thereof) and from above by the distance of 
the classical cyclic code. Since GV bound always pro- 
duces asymptotically good codes, the parameters of our 
quantum codes will be mostly limited (at least, for long 
blocks) by their binary counterparts. 

We begin with analyzing the condition ( p2] ) for a cyclic 
CWS code. For a vector e 6 FJ to be in C±, this con- 
dition can be rewritten in terms of the corresponding 
polynomials, 

p(x) [u(x) + r(x)v(x)] = modi" - 1 , (36) 

where the coefficients of the (reversed for notational con- 
venience) polynomial e(x n ~ 1 ) = u(x) + uiv(x) arc given 
by the components of the vector e e FJ. Since binary 
p(x) divides x n — 1, we can rewrite this in terms of the 
corresponding generator polynomial q(x) = (x n — l)/p(x) 
for the binary code C, 



(37) 



u(x) + r(x)v(x) = modq(x). 



Now, if v(x) is mutually prime with q(x), Eq. ( |37j ) can 
be just solved for r(x). In this case the answer is unique 
[modq(a;)]. On the other hand, multiple solutions for 



r(x) are possible when gcd[v(x),q(x)] ^ 1. In this work, 
we avoid the complications caused in the latter case [[33) 
and only consider irreducible polynomials q(x). 

Overall, for any irreducible q(x) and any e with v/0 
and wgt e < d(C), Eq. ( |37j ) has a unique solution for r(x) 
such that degr(x) < degq(x) = n — k. Respectively, 
there is no more than one additive quantum code such 
that e 6 C±. Generally, only some of thus obtained r(x) 
correspond to self-orthogonal codes, see Eq. (|33|). 

Below we complete the greedy argument by counting 
the polynomials r(x) corresponding to self-orthogonal 
codes, Eq. (|33|). We consider separately the case 
when the irreducible polynomial q(x) is palindromic [see 
Eq. (|35])] in Lemma |t] below, and when it is not in 

Lemma 5. Consider a cyclic binary code C[n, k, de] with 
the generator polynomial q(x) which is both irreducible 
and non-palindromic, x degq ^q(x~ 1 ) ^ q(x). Then there 
exists a single- generator additive cyclic code [[n, k, d]] 
with distance 



d = mm(dc-, d, 



GV) 



restricted by both the distance dc of the binary code and 
the following variant of GV bound 



d-l 



d GV = maxd:£(3 s -3)I*Z^V 



s=l 



(38) 



Proof. The non-palindromic generator polynomial q{x) is 
one of the factors in x n — 1 which also contains its recipro- 
cal, x desq ^q(x~ 1 ). This implies that the corresponding 
parity-check polynomial p(x) satisfies Eq. d35|). Further, 
since q(x) is irreducible, the solution r(x) of Eq. ( |37| ) 
is unique assuming v(x) ^ and degr(x) < n — k, 
which gives the exponential term in the r.h.s. of Eq. (38). 



Eq. (38) improves on the standard GV inequality (27) by 
discarding a few sets of vectors. The first set are vectors 
with u(x) = 0, which implies r(x) = modg(a;). The 
second set are vectors with u(x) = v(x), which all give 
r(x) = 1 modq(x). The third set are non-zero vectors 
with v(x) = 0, which can never be in Cj_ corresponding 
to the generator (^o|) . Finally, note that any error vector 
of weight s produces at least n/ gcd(s,n) different cyclic 
shifts. All of these cyclic shifts give the same polynomial 
r{x) and can be discounted. The condition d < dc comes 
from Lemma |l|. □ 

Note. The lhs of bound ( |38|) limits the number of cyclic 
classes for 4-ary vectors e of weight s < d — I. Most of 
these vectors have the maximum period n. Therefore, it 
can be proven that the term gcd(s,n)/n in (38) can be 
replaced with the smaller term that rapidly tends to 1/n 
for large n. In turn, boun d (|38| ) adds about log 2 n infor- 
mation qubits to bound (|27|)but tends to the standard 
quantum GV bound ( |23| ) as n — > oo. 

In the following example, this bound coincides with 
inequality d < dc, which uniquely sets code distance d. 
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Example 8. The family of the binary codes with the 
parameters [n = 2 4h + 2 3h - 2 h - 1, k = n - 6h,3], 
h = 1,2,..., constructed in Ref. ^5^, has irreducible 
non-palindromic generator polynomials as required in 
Lemma |^. For dc = 3 the sum in Eq. ( jgq j has just 
one term at s = 2; for n odd gcd(n, 2) = 1. Explicit cal- 
culation confirms that the parameters of these codes sat- 
isfy inequality (f?is[), which proves the existence of single- 
generator cyclic quantum codes with exactly the same 
parameters, [[n = 2 4/l + 2 3h - 2 h - 1, k = n - 6h,3]}, 
but not necessarily cyclic CWS codes. The smallest 
of these codes, [[21,15,3]], corresponds to a polynomial 
q(x) = 1 + x + x 2 + x 4 + x 6 (unique up to a reversal) and 
can be obtained from an order-21 circulant graph corre- 
sponding to r(x) = x+x 4 +x 17 +x 20 . This particular com- 
bination of parameters gives the best existing codeWS^. 



Example 9. According to the BCH bound J56| /, a cyclic 
code has distance d>r+l(r + lis the "designed" 
distance) if the corresponding generator polynomial q(x) 
has r consecutive roots, e.g., a, a 2 , . . . , a r , where a is 
the primitive n th root of unity. A polynomial m a (x) 
which has root a, necessarily has s distinct roots a 23 for 
all j = 0, s — 1 if s is the smallest number such that 
2 s = 1 mod n. We say that the code has zeros a 1 , where 
exponents i form the set I — {2 J (modn), j = 0, s — 1}. 
The code generated by m a {x) has designed distance 5 if 
3 £ J or, equivalently, if 2 s = 3mod?i for some s. The 
polynomial m a (x) is non-palindromic if — 1 ^ I. We can 
further obtain codes with irreducible non-palindromic 
generators and designed distances 7, 9, etc., by imposing 
additional conditions, e.g., 5 E 1 , 7 E I , etc. The values 
of n, for which this is possible form an infinite set, 
{23 5 , 47 5 , 71 7 , 95 5 , 115 6 , 143 B> 167 5j 191 7 , 235 5 , 239 7 , . . .}, 
where the subscripts indicate the designed distances. 
In fact, the first three codes represent the well known 
quadratic-residue codes with the higher distances (exceed- 
ing the BCH bound) equal to 7, 11, and 11, respectively. 
GV bound proves the existence of additive quantum CWS 
codes with the parameters [[23,12,4]], [[47, 24, d > 6]], 
[[71,36,d > 7]], [[95,59,5]], [[115,71,5]], [[143,83,11]], 
etc. The first three codes have the parameters as good as 
any known codes with such n and k. 

Now let us consider the case of a palindromic polyno- 
mial q(x). First, we prove 

Lemma 6. Consider a binary code C generated by a 
palindromic polynomial q(x) such that q(l) = 1. Then 
any quantum code $3(\) which satisfies self-orthogonality 
condition (ffi|) is equivalent to a cyclic CWS code with a 
symmetric polynomial r{x), see Eq. (ffi[). 

Proof. The corresponding check polynomial p(x) is sym- 
metric, thus the condition (B3[) can be rewritten as 



r(x) + r{x 



,n-l- 



mod q(x). 



(39) 



The condition q(l) = 1 guarantees that the palindromic 
polynomial q(x) has odd weight and even degree 2m, in 



which case the "central" monomial x m has non-zero co- 
efficient q m = 1. Given the block length n, let us choose 
an equivalent code [see Eq. (|32|)] with r(x) such that the 
coefficients 



Tm+l — r m+ 2 — 



= 0. 



(40) 



The coefficients of the polynomial in the l.h.s. of Eq. ( J39| ) 
satisfy the same condition, except for the term x n ~ m , 
which has coefficient r m . The coefficients are arranged 
in such a way that the l.h.s. of Eq. (||^) can only equal 
zero or x n ~ m q(x) mod x n — 1. However, the latter pos- 
sibility can be excluded by comparing the corresponding 
coefficients of the free term x . The only remaining case 
corresponds to a symmetric r{x). □ 

It is now clear that for a palindromic irreducible gen- 
erator polynomial q(x) ^ 1 + x, one should reduce the 
count in the r.h.s. of Eq. (^8|) by replacing 2 n ~ k with 
2("~ fc )/ 2 j the number of symmetric polynomials that sat- 
isfy Eq. (40). This gives the following version of GV 
bound, 



^gv 



axd : ]T (3 s - 3) gCd( J'" } Q < 2^' 2 2. 

(41) 

While the resulting estimate is much weaker than the 
GV bound (|3|), it still gives asymptotically good codes. 
A better (especially, for small d) bound is given in the 
following 

Lemma 7. Consider a cyclic binary code C[n, k,dc] 
with dc > 3 and the generator polynomial q(x) which 
is both palindromic and irreducible. Then there ex- 
ists a cyclic CWS code [[n, k, d}] with the distance d = 
min(dc, [_<^Gv/2j) 7 where day = maxcZ: 

y (3 r./*i _ 3) i£dM (>/2j\ < 2(n _ fc)/a _ 2 (42) 



Proof. The restriction on the distance guarantees that 
q(x) 7^ 1 + x, and therefore q(x) satisfies the conditions 
of Lemma ^; in particular, n — k is even. The inequality 
just corresponds to symmetric polynomials r(x) and the 
errors that are also symmetric, e(x n ~ 1 ) = e(x) modi" — 
1. The statement of the Lemma follows from the fact 
that for any general error e(x), there is a symmetric error 
(x) = e(x) + e(x n ~ 1 ) modi 11 — 1 whose weight is even 

□ 



-'sym 



and is limited by wgte sym (a;) < 2wgte(x) 



Example 10. Among classical codes, the largest distance 
is obtained for the repetition codes, with the parameters 
C = [n, 1, n] . The parity-check polynomial is p(x) = x+l; 
the generator polynomial q(x) = 1 + x + . . . + a;" -1 is 
irreducible and palindromic for n = 2, and for all n > 2 
that satisfy the condition ord m (2) = m— 1, where ord m (q) 
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is the multiplicative order of q modulo m. 
the following n < 100: 



This includes 



{3, 5, 11, 13, 19, 29, 37, 53, 59, 61, 67, 83, . . .}. (43) 

Lemma [? shows that for n from the set (^^), additive 
cyclic CWS codes with parameters [[n, 1, dny(n, 1)/2J]] 
exist, where dcvin, 1) is obtained from Eq. ( Ua j with k = 
1. Asymptotically, at large n, this corresponds to cyclic 
codes with the relative distance given by half of that given 
by Eq. (|§j. 

Example 11. (Cyclic analogs of the toric code) In 
the setting of the previous example, cyclic CWS codes 
[[5,1,3]], [[13,1,5]], [[25,1,7]], [[41,1,9]] with p(x) = 
1 + x were obtained numerically. The corresponding 
graph-state generators are ZXZ for d = 3, ZZXZZ for 
d = 5, etc. We obtain a family of cyclic codes with the 
weight-A generators, S3 = ZXXZ , S$ = ZIXXIZ, etc. 
Codes with generators S3 = ZXIXZ , S5 = ZXIIIXZ , 
S 7 = ZXIIIIIXZ, and S 9 = ZXIIIIIIIXZ have the 
same parameters (the corresponding graphs are some- 
what more complicated). The latter family can be gen- 
eralized to codes with n = t 2 + (t + l) 2 , k = 1, d = 
2t + 1, t = 1,2, .. .; the corresponding stabilizer gener- 
ators $24+1 having It — 1 identity operators separating 
ZX and XZ. These cyclic codes are related to gener- 
alized toric codes t ~5o\ ~5%; the square-lattice qubit layout 
preserving the circulant symmetry is illustrated in Fig. Q 
fort =1,2. 
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FIG. 3. (Color online) Left: Correspondence between the 
cyclic code [[5,1,3]] with generators ZXIXZ, and a gen- 
eralized toric code on square lattice. Only qubits within 
the dashed square are participating in the code; periodic 
boundary conditions are assumed. The vertical dashed line 
indicates a topologically nontrivial chain of errors which 
limits the distance of the code: is equivalent to 

ZxZ^Z^Z^Z^,. Right: same for the code [[13,1,5]] with the 
generator ZXIIIXZ. 



Example 12. (CWS codes from k copies of the repeti- 
tion code ) Let us take the binary code C to be formed by 
k copies of the repetition code with the distance c?2 = m. 
Then the block size is n = km, and the check poly- 
nomial is p(x) = x k — 1. The generator polynomial 
q(x) = 1 + x k + . . . + a^" 1-1 ' is always palindromic; 
it is also irreducible if m belongs to the set (45), while 



k = m s , s = 0, 1, 2, . . .. For sufficiently large n, Lemma^ 
gives asymptotically good codes with kd oc n 2 . Since 
these parameters cannot exceed those of the binary code 
C which correspond to kd = n ( thus 6 = 1/k), for these 
values of m and k > 10, there exist cyclic CWS codes 
with the parameters of the corresponding binary code, 
[[n = m s+ ,m s ,d = m]]. This prediction is readily veri- 
fied empirically, see Table [|. Note that, as in the Example 
[71| , many of these codes have stabilizer generators with 
small weight. 



m 


n 


k 


d 


Si 


3 


6 


2 


2 


ZXZ 




9 


3 


3 


ZXZ 




12 


3 


3 


ZXZ 


5 


5 


1 


3 


ZXZ 




10 


2 


3 


ZXZ 




15 


3 


5 


ZIZIIXIIZIZ 




20 


4 


5 


ZIZZXZZIZ 




25 


5 


5 


ZIZZXZZIZ 


7 


7 


1 


3 


ZXZ 




14 


2 


5 


ZZIXIZZ 




21 


3 


6 


ZIZZXZZIZ 




28 


4 


7 


ZIZZIXIZZIZ 



TABLE I. Families of the cyclic codes obtained numerically 
from k copies of the classical repetition code, p(x) = x m — 1, 
corresponding to m = 3, 5, 7. The expected distance satura- 
tion, d = m, is reached already at k < 4, even for k and m 
different from those in Example [l^. The operator strings in 
the last column are representative graph-state generators (the 
remaining generators are obtained by cyclic shifts). 



VII. CONCLUSIONS 

In this paper, we analyze how the general CWS frame- 
work can facilitate the search of the additive quantum 
codes with reasonably good parameters. Unlike complete 
optimization of CWS codes |2(J, which involves all non- 
isomorphic LC-inequivalent order-n graphs and all binary 
codes of length n, here one can independently pick a suit- 
able graph Q, and then search for a linear binary code C 
that can correct the error patterns induced by Q, see 
Eq. pi). 



The choice of the graph is discussed in Sec. [II. In the 
simplest case of pure codes, one has to pick a graph with 
a sufficiently large graph-state distance d'(Q). Assuming 
that a regular graph is being sought in this design, we 
consider graphs with minimal vertex degrees d'(G) — 1 or 
more. 

After the graph is chosen, the second step involves the 
search of an appropriate binary code. Here in Sec. [v] we 
prove the existence of the binary codes that give good 
quantum CWS codes. The corresponding lower bound on 
the distance is given by the quantum Gilbert- Varshamov 
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bound (p7|). Note that while this bound is proved for a 
given graph, it is the same bound that holds for a generic 
stabilizer code. 

Our results show that by restricting the graph Q of a 
CWS code to regular lattices, one can lower the com- 
plexity of the code search and still obtain codes with rel- 
atively good parameters. On the other hand, the graph 
structure could be mapped directly to a physical qubit 
layout. Therefore, such codes can simplify both the hard- 
ware design and the error-correcting procedures, which 
will easily admit the property of translational invariance. 

An unexpected byproduct of this work is the discov- 
ery of a previously unexplor ed fam ily of single-generator 
quantum cyclic codes (Sec. |VI A ). These codes are rel- 
atively easy to construct, and they are plentiful. We 
construct (or prove the existence) of several simple fami- 



lies of such codes that have unbounded distances. These 
include cyclic CWS codes with weight-4 stabilizer gen- 
erators, which turned out to be toric codes in disguise 
(Example |ll| ) , as well as a code family with the parame- 
ters of generalized repetition codes, [[kd, k, d]] (Example 
|l2| ). The main advantage of these families is a simple 
structure of their stabilizers. 
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